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We have found the exact (factorized) ground state of a general class of ferrimagnets in the presence 
of a magnetic field which includes the frustrated, anisotropic and long range interactions for arbitrary 
dimensional space. In particular cases, our model represents the bond-alternating, ferromagnet- 
antiferromagnet and also homogenous spin s model. The factorized ground state is a product of 
single particle kets on a bipartite lattice composed of two different spins (p, a) which is characterized 
by two angles, a bi-angle state. The spin waves analysis around the exact ground state show two 
branch of excitations which is the origin of two dynamics of the model. The signature of these 
dynamics is addressed as a peak and a broaden bump in the specific heat. 

PACS numbers: 75.10.Jm, 03.67.-a, 64.70.Tg 



Introduction.- Spin models are the building blocks of 
the theory of quantum magnetism and strongly corre- 
lated electron systems. In addition, they have been con- 
sidered as an effective model to describe the behavior 
of a system in several disciplines. Recently, the im- 
plementation of quantum notions in quantum devices 
has attracted much attentions both in research labs 
and demanding applications like nanotechnology, quan- 
tum computation]]! and particularly optical lattices 0. 
Quantum spin models arc prototype realization of many 
relevant properties of quantum implementation in such 
devices. Therefore, different aspects of a quantum 
phase is of utmost importance for scientists and engi- 
neers. Quantum phases are characterized by the ground 
state (GS) properties of the corresponding many body 
system. [U 

Except of a few particular cases such as ID bond al- 
ternating Heisenberg anisotropic Heisenberg model 
(XYZ), XXZ in a longitudinal magnetic field and Ising 
model in transverse field which are exactly solvable 
the GS of a general spin model is not known. However, 
at some special values of the model parameters the quan- 
tum correlations are vanishing and the GS is a product 
of single particle states. The factorized state (FS) man- 
ifests zero entanglement which is necessary to be iden- 
tified for reliable manipulating of quantum computing. 
A FS (unentanglcd) which is associated with an entan- 
glement phase transition can be also a quantum critical 
point in certain condition which is discussed in this arti- 
cle. This information is also attractive for the study of 
quantum phase transitions. Moreover, finding an exact 
ground state (as a FS) even at particular values of the 
parameter space of a many body spin model leads to the 
identification of that phase in addition to more knowledge 
about the properties of the model close to the factorized 
point via implementing an approximate method. 



In a seminal work, Kurmann, Thomas and Miiller [6( 
identified the factorized state of a homogeneous spin-s 
XYZ chain at a magnetic field of arbitrary direction. 
Factorized GS has been also observed in the two di- 
mensional lattice through Quantum Monte Carlo sim- 
ulation in terms of entanglement estimators. Recently, 
Giampaolo, Adesso and Illuminati Q have introduced a 
general analytic approach to find the factorized ground 
states in a homogenous translational invariant spin-,s 
quantum spin model for arbitrary long range interaction 
and any dimensional space. Their study is based on the 
single-spin unitary operation and the factorized point is 
determined at the position where the associated entan- 
glement excitation energy becomes zero. More elaborate 
explanations came up with some generalizations quite 
recently. The factorized GS of the dimerized XYZ spin 
chain in a transverse magnetic field has been investigated 
and reported that the factorized point in the parameter 
space of the Hamiltonian corresponds to an accidental 
ground state degeneracy, (lfjj However, in this article we 
will present (i) the FS of an inhomogencous (ferrimag- 
netic) spin model which is composed of two spins (p, a) in 
the presence of a magnetic field on a bipartite lattice with 
arbitrary long range interaction and dimensional space, 
(ii) the Hamiltonian is not necessarily translational in- 
variant and (iii) the exchange couplings can be compet- 
ing antifcrromagnctic and ferromagnetic arbitrarily be- 
tween different sublattices to build many practical mod- 
els such as frustrated, dimerized and tetramcrized ma- 
terials. Moreover, our results recover the previous ones 
for a = p and a particular configuration of the couplings. 
0, 0, B El ln addition, we will address on the existence 
of two energy scales which lead to a surprising dynam- 
ics of the model close to the factorizing point and its 
fingerprint as a double peak in the specific heat versus 
temperature. As an enclosure, the results have been ap- 
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plied to the ID ferrimagnctic XXZ (p, a) spin chain in 
the presence of a transverse ma gne tic field which is real- 
ized as a bimetallic substance, [ljj We will also address 
the cases where the factorizing field coincides the critical 
point. 

Factorized state.- Let us consider a two sites model 
which is composed of two spins a = and p = 1 with 
the following Hamiltonian 



Tjf jX X -x 

rl = J a p 



J v <7 v p v + J z a z p z + h'(a z 



(1) 



where J^,p = x,y,z are the exchange couplings in dif- 
ferent directions and h! is proportional to the magnetic 
field. We are looking for a factorized state which is sat- 
isfied by H'\a)\p) = e\a)\p), in which |er) and \p) arc 
the single particle states. It is appropriate to choose \a) 
and \p) to be the eigenstates of a • n' and p ■ n" with 
eigenvalues +| and +1; respectively, where n'(9,ip) and 
fi"(f3,a) are unit vectors in Bloch sphere. The solution 
of H' | cr) | p) = e\a)\p) gives the factorized state at h' = h'f 
and its corresponding energy (e) (l2l | . Moreover, we found 
that the angles 9 and (3 are fixed by the couplings (J^ 7 h'; 
see Eq.@) while a and ip are given by one of these choices 
(I) a = 0, cp = 0, (II) a = 0, <p = tt, (III) a = § , ip = -f , 
(IV) a — ^ , ip = ^ . The spins are located in the xz-planc 
for choices I and II while they have projections only in 
the yz-plane for III and IV. Without loss of generality we 
can assume the spins are located in the xz-plane. In fact, 
the spins of yz-plane will fall to xz-planc by interchange 
of J x <-> J y . Moreover, the coordinates (0,<p = 0) and 
(— 0, tp = tt) are representing the same direction, there- 
fore case (I) a = 0, ip = is able to describe all possibil- 
ities. 

The two spin model (cr = h, p = 1) is now general- 
ized to arbitrary (a, p) spins, [y] To find the factorized 
state of a general two site ferrimagnet we consider a ro- 
tation on a and p spins such that ~a and p point in 
(8, tp = 0) and (/3, a = 0) directions, respectively, The 
rotation operator is D = D a (0,8,0)D p (0, (3,0) where 
DP(O,0,O) = D(a = 0,/3,7 = 0) = D z {a)D y {l3)D z (~f) 
is defined in terms of Eulcr angles and a similar expres- 
sion is considered for D a (0, 6>,0). Then, We impose the 
condition to have a factorized (fully polarized) eigenstate 
for this Hamiltonian which fixes the following relations 
for the model parameters 



cos ( 



cos (3 



h', = 



h'fjy + J X (J Z ~ - jy )pa + h' f J z (jyp + J x a) 
h!fj x + jy(J z2 - J x2 )pa + h' f J z (J x p + jy<j) ' 

h'fjy + j x (.j z2 - jy 2 ) P a + h' f j z (jy<j + j x p ) 

hfJ x + Jy(J z2 - J x2 )pa + h' f J z (J x a + Jvp) ' 
-{2J x JVpa+ (p 2 + a 2 )J z2 +CJ Z ), 



C = yj4pa(pJ x + aJ y ){aJ x + pjv) + (p 2 - a 2 ) 2 J z \ 

j x jy 



(2) 



Therefore, for arbitrary (a, p) and at the above value for 
h' = h's we have a fully polarized eigenstate which is a 
factorized state. The ordering of this state is defined by 
two angels (9, j3) which show the orientations of ( a , p ), 
respectively. 

Now, we intend to find the condition for having a fac- 
torized state for a ferrimagnctic lattice in a magnetic 
field. We will then show that the factorized state is the 
ground state of lattice implementing a spin wave theory 
to study the quantum fluctuations. We consider a general 
Hamiltonian of ferrimagncts on a bipartite lattice where 
sublattice {A a ) contains a spins and the other sublat- 
tice (Bp) includes p spins. The interaction can be long 
ranged between different sublattices but no interaction 
in the same sublattice. The ferrimagnctic Hamiltonian 
for such case can be written as 



H = 



E 



^CiCi+r(t/ r Cj Pi+r + a i Pi+r) 



JZ z z 
'r °i Pi+r 



(3) 



where i = 12,^3) and r = (r±,r2,Ts) are repre- 
senting the three dimensional index on the lattice and 
Ci,Q+r = ±1 which realize both ferromagnetic (F) and 
antifcrromagnetic (AF) exchange interactions. A remark 
is in order here, the Hamiltonian in Eq.([3]) is a sum of 
two sites Hamiltonian defined in Eq.([T]) where the two 
spins can be far from each other. However, the inter- 
action between each couple of ((Tj, Pi+ r ) can depend on 
distance (r) with different strength and also be F or AF 
arbitrarily defined by Q, Ci+r- A factorized eigenstate for 
the Hamiltonian of Eq. Q can be written as 



\FS) 



(4) 



where |c^) and are the eigenstates of <?, •h' i and pj - h" 
with largest eigenvalue where f)! i and fi" are unit vectors 

pointing in (Q0,(p = 0) and (Qf3,a = 0), respectively. 
However, the factorized state d-FS')) is an eigenstate of 
the Hamiltonian if the angle Q9 (Cj(3) be consistent with 
all pair of interactions originating from at (/?,) on sublat- 
tices A a (Bp). According to Eq.© the former condition 
is satisfied if the interaction between each pair (at , Pi+ r ) 
is the same for all directions while depending on distance 
(r), i.e = A(r)J^,/i = x,y,z,X(r) > 0. Under these 
constraints the factorized state (Eq.Q) is an eigenstate 
of H with the characteristic angles (9, (3) defined in Eq. j2]) 
and the factorizing field is 



(5) 



where N r is the number of spins on each sublattice. To 
show that \FS) is the ground state of H at hj we first 
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implement a rotation on the Hamiltonian. The rotated 
Hamiltonian (H) is the result of rotations on all lattice 
points of H, 

H = D^HD ,D= (g) Df(0, Q9, 0)Z>?(0, Q(3, 0). 
ieA a ,jeB p 

(6) 

In the next step the rotated Hamiltonian is bosonized 
using the Holstein-Primakoff (HP) transformation, erf = 



'2(7 - ojoj a,, of = a-a]a i ,Pj' = yj2p — bjbj bj, pj = 

p — b^bj, where ai(a\) and bj(bj) arc two types anni- 
hilation (creation) boson operators. Using HP trans- 
formations the anisotropic ferrimagnetic spin model is 
mapped to an interacting system of bosons. The Hamil- 
tonian in the momentum (k) space and in the linear 
spin wave theory is diagonalized via the rotation, Xk = 
a k cosf] k - e lS h sin rj k ;i' k = e l6 b k cosij k + a k sinr) kl and 
a shift at k = where 6 is defined by J2r e ~ lkr Jr = 
| e ~ lkr Jr\ elS ■ The diagonalized Hamiltonian is 

H = E gs + (k)xlxk + uj + (k)^kj , (7) 

where E gs is the ground state energy [T^j and o; ± (/c) > 
are normal modes parallel and perpendicular to the field 
direction. 



,t. 



: {k) = D+ ± 



D- + V^tan(27 7fe )| £ r e^*"' J* 



tan(2ry fe ) 



v /l + tan 2 (2r ?fc ) 
D- 



(8) 



in which 



D ± = h-£±l) + h f (fc OS 9±£- cos/3) 
2r z p u 2p 2<t 

r 1 ^ hf — h 

- -^ r (<J±p) + ^—(cos(3±cos9), 

= J"^A(r). 



(9) 



The bosons number (<ra) is proportional to (hf — h) 
which states that at h — hf the bosons number in the 
ground state is exactly zero [1^ . i.e (ojcit) \h=h f = = 
{b\bi)\h=h r It is also true for (a\bi)\ h = hj = 0. Thus, 
the spin wave theory is exact in all orders at the factoriz- 
ing field [h — hf) and \FS) is the corresponding ground 
state of H . 

To visualize the configuration of a factorized ground 
state of a general frustrated model we have plotted an 
example in Fig. (TT]) with the assumption J^,J^ > and 
J,f < where Q and Q +r define the sign of interac- 
tions. The solid-lines represent antiferromagnctic inter- 
action and the dash-dotted ones are the ferromagnetic 
counterparts. As shown in Fig.(TT]) the interactions can 




FIG. 1: The configuration of a factorized state on a one- 
dimensional lattice for arbitrary frustrated couplings. Solid 
lines (dash-dotted) represent antiferromagnetic (ferromag- 
netic) couplings which are defined by QQ+r as depicted by 
± on each site. Each color belongs to equal distance interac- 
tion (same r). 



be frustrated and long ranged without a translational in- 
variancc. However, the factorized state is defined by two 
angles (6,(3) while each a (p) spin is directed in 9 ((3) 
or —9 (~(3) directions. We call this a bi- angle ordering. 
In a special case the bi-angle ordered state can configure 
a ferromagnet (6 = (3) or antiferromagnet (9 = it — (3) 
factorized state. 

Discussions.- In a spin model when the magnetic field 
is strong enough all spins will align in the direction of the 
magnetic field which characterizes the saturated phase as 
far as h > h s . In our notation, the saturated phase ap- 
pears when all a (p) spins get 9 — tt ((3 = tt). Thus, the 
saturating field (h s ) is a factorizing one when J x = J v . 
In case of J x ^ J y the saturation can only appear at in- 
finite value of the magnetic field while a finite factorizing 
point (hf) still exists. For J x = J y , the lower excitation 
band becomes gapless (uj~(k = 0) = 0) which confirms 
that the factorizing point (hf = h s ) is the critical point 
which separates the non-saturated phase (h < h s ) from 
the saturated one (h > h s ). It is worth to mention that 
at J x — J v the rotational symmetry around the magnetic 
field is restored where the quantum fluctuations around 
the field axis are suppressed. In the absence of rotational 
symmetry the saturation in the field direction does not 
happen at a finite value of magnetic field while the model 
approaches saturation asymptotically at infinity. 

A general feature of our result is that it can simply 
recover the previous study of homogenous systems by 
replacing a = p = s. In that case the restriction of 
bipartite lattice is promoted to arbitrary lattice and the 
interaction between any pair of spins can exist. However, 
our Hamiltonian is not restricted to the translational in- 
variant symmetry or bond-alternating ones which is wit- 
nessed by the example given in Fig.{T]). This can also 
be generalized to any dimension. We claim that the gen- 
eral Hamiltonian which can possess a nontrivial factor- 
ized ground state should be of the form Eq.([3]) with the 
restriction = X(r) J M , p = x, y, z. 

In the absence of magnetic field, the spin-(/9, a) ferri- 
magnet turns into the spin-cr antiferromagnet in the limit 
p — > cr, whereas it looks like the spin-p ferromagnet in the 
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other extreme limit £ — » oo. In this sense, the difference 
p — a can be regarded as the ferromagnetic contribution. 
The analysis of our result shows that the value of hf is a 
decreasing function of £ converging to the ferromagnct 
feature as this ratio becomes large. More investigations 
on different applications of our approach in several mod- 
els are in progress. [l2| 

Let us now be more concrete by concentrating on the 
one-dimensional nearest neighbor XXZ ferrimagnet in 
the presence of transverse magnetic field. Suppose that 
jx _ jz _ j anc j jy _ j ^ wnere A represents the easy 

axis anisotropy. At zero magnetic field the quantum fluc- 
tuations are large and the ground state of the model is 
strongly entangled. Upon adding the transverse magnetic 
field the U(l) symmetry of the XXZ model is lost and 
the entanglement of the GS is decreased. In the mapped 
bosonic system the magnetic field is served as a chemical 
potential, thus the number of bosons ((at a) and (tfb)) is 
dependent on the magnetic field. An enchantment of the 
magnetic field causes deducing of the bosons' number and 
the quantum correlations decrease. At factorizing field 
h = hf, the number of bosons is zero and the quantum 
fluctuations become completely uncorrelated. Moreover, 
our calculations show that the factorizing field in fer- 
rimagnetic model depends on the anisotropy parameter 
(A) similar to a homogeneous antiferromagnetic Hciscn- 
berg model. For A = (fcrrimagnetic XY model) the 
in-plane magnetic field completely breaks the rotational 
symmetry Increasing A suppresses the effect of the field 
and try to evoke the rotational symmetry to the sys- 
tem. Thus, by increasing A the factorizing Neel field 
approaches to the saturation field. At A = 1, the rota- 
tional symmetry is completely repayed and hf is exactly 
lied on h s . 

A benefit of identification a factorized state is that we 
can work out an approximate method around the factor- 
ized point to get some information on the properties of 
that phase. This help us to calculate the magnetic prop- 
erties of the ferrimagnetic XXZ model in the presence of 
a transverse magnetic field. We have implemented the 
linear spin wave theory around h = hf for a = ^ and 
p = 1. Our results for magnetization (M x , M y ) and stag- 
gered magnetization (SM X , SM y ) in both x and y direc- 
tions arc plotted in Fig. (J3J) where the magnetic field is in 
x direction. The approach can be extended to derive the 
thermodynamic of this model which will be done in near 
future. 

It is also worth to mention that our results are ap- 
plicable to the homogenous XXZ Heiscnbcrg spin-1/2 
chains in the presence of a transverse magnetic field 
(h x ). This model has been studied intensively in the 
literature. [l3, 14, ID, Ii| On the onset of transverse mag- 
netic field a perpendicular anti-ferromagnetic order is sta- 
bilized by promoting a spin-flop phase. At the classi- 
cal (factorizing) field (hf < h c ) in the spin-flop phase 
the ground state is factorized to single spin states and 



1.4 


: - Mj 






1.2 


: - - M > 






1.0 


: " - SM X 


\ / 




0.8 


; . SM V 






0.6 








0.4 








0.2 
0.0 




s __ 





FIG. 2: The magnetization and staggered magnetization of an 
anisotropic ferrimagnetic (a — 1/2, p = 1) spin chain versus 
transverse field and for A = 0.25. 



the staggered magnetization along the y direction has a 
large value. Very close to the critical field (h c ) the anti- 
ferromagnetic order becomes unstable and the staggered 
magnetization falls sharply to vanish at the critical point. 
For h > h c the spins are almost aligned in field (x) direc- 
tion and the factorized state as a fully polarized phase 
will be appeared for h — > oo. The excitation energies 
around the factorizing point have the following form: 



w ± (fc) = (l + A)A + A(±cos( ^ ) 



!)■ 



(10) 



Thus, we have two branches of magnon energies as two 
scales which impose two dynamics in the system. The 
most interesting feature is that around the factorizing 
field both scales show up. These dynamics correspond 
to the coexistence of two different features of the model. 
In other words, by increasing the transverse magnetic 
field, spins try to align in the x direction which is the 
ferromagnetic feature of the system. In the intermediate 
values of h, we have both anti-ferromagnetic and ferro- 
magnetic behaviors since the anti-ferromagnetic order in 
y direction has already been stabilized. The finger print 
of these features appear in the thermodynamic functions 
such as s pec ific heat and internal energy. As it is seen 
from Ref. [16| the second feature can be seen as a shoul- 
der at the right side of specific heat curve. By further 
increasing of h, the ferromagnetic behavior is seen as a 
broaden peak in the curve. This point is almost near the 
classical field where the ground state of the system has 
been factorized. 
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Excellence in Complex Systems and Condensed Matter 
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